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1. Introduction 

We performed simulations for 2 + 1 light flavors Lattice QCD with the fixed-point Dirac oper- 
ator. The (parametrized) fixed-point operator we use approximatively satisfies the Ginsparg-Wilson 
relation [Hp: 

Dy 5 + Y 5 D = DY 5 2RD , (1.1) 

and is limited to a hypercube [Q]. (Here R is a non-trivial local matrix.) The quark mass is intro- 
duced by 

D{m)=D + m[^- l -D ] j. (1.2) 

The fixed-point action gave very promising results in the quenched approximation, in particular 
good scaling even at a = 0.15fm ^J. We tuned the coupling in our full QCD simulation to be 
close to a = 0. 15 fm [||], and measurement of the Sommer parameter ro gave a result in agreement 
with this value. We simulate two lattices, of size 8 3 x 24 and 12 3 x 24 at this coupling. The results 
presented here refer to the smaller lattice. 

2. Autocorrelation times 

We estimated the autocorrelation time of plaquettes built from smeared links and that of the 
pseudoscalar-pseudoscalar correlator. We are using a global update [§], the configurations are 
separated by « 0.7 standard Metropolis sweeps. 

Fig. [I] shows the MC history for the plaquette for a longer run and the autocorrelation function 
extracted from it. This gives an autocorrelation time To ~ 40. 

In Fig. |2| we plot the autocorrelation function for the zero-momentum (P(0)P(t)) correlator at 
different time separations t. The autocorrelation time for this quantity is estimated to be somewhat 
higher. 
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Figure 2: The autocorrelation function for the (P(Q)P(t)) correlator at different t values. 



3. The delta regime 

The epsilon regime is defined by: 



m PS L, m PS L 4 < 1 



(3.1) 



and the box size should be large enough, 2nFL ^ 1, for the validity of chiral perturbation theory. 
Here mps is the pseudoscalar mass in infinite volume. 

In the delta regime [§] one considers a lattice under similar conditions but elongated in the 
time direction so that the low-lying spectrum in the corresponding spatial volume can be measured. 
In this regime, for two massless quarks, the low lying spectrum is given by an 0(4) rotator having 
in leading order in 1/L a mass gap 

3 

(3.2) 



2F 2 L? 



Note, however, that for our box size of L = 1 .2 fm one expects non-negligible corrections to this 
leading order result |0]. 



4. Low-mode averaging 



We used the method of low-mode averaging |§ |, [K]] to decrease the error in the zero- 
momentum meson correlators 



C(t;y) = (^{n75G(x,y)h5r 2 G(x,y)) 



(4.1) 



t=X4-y 4 



where F\ , T2 are arbitrary Clifford algebra matrices. As part of our updating algorithm we calculate 
the low-lying eigenvalues and their corresponding eigenvectors ^ which allows us to perform such 
averages. For the 8 3 x 24 lattice we stored 48 eigenvectors with the smallest eigenvalues. We split 



up the correlator in Eq. (1.1) into three parts: 

C(t;y) = Cu{t)+C lh {t;y Q ) + C hh {t;y) , 



(4.2) 
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# confs 

Figure 3: The zero-momentum (P(O)P(f = 5)} correlator measured at fixed source position (Co); by using 
the eigenvectors to average the source position over the whole lattice in the low-low part (Ci); by averaging 
in the source position over a time-slice in the low-high part as well (C2). 



where / refers to the low-mode part Gi(x,y) of the quark propagator (which is given by the calcu- 
lated eigenvalues and eigenvectors), and h to the rest of the propagator, G(x,y) =Gi(x,y)+Gh(x,y). 
We analyze the effects of low-mode averaging for Qi and Cy, separately. 

On Fig.|| we compare the zero momentum correlator (P(O)P(t)) at time separation t = 5 for 
the different amounts of averaging. Co is the correlator without any averaging over the source 
position, C\ is when the 11 part is averaged over the whole volume, C2 is when in addition the Ih 
part is also averaged over a time slice. While the // averaging is very effective, the Ih averaging 
barely shows any suppression of the fluctuations. Since the latter one is also more expensive, one 
can safely ignore this option and average only in the // part. 



5. Results 



5.1 The PCAC mass 

Because we are using the parametrized fixed-point action which solves the Ginsparg-Wilson- 
equation only approximately we get an additive mass renormalization for our PCAC masses. The 
size of this mass shift is an indicator for the quality of our approximation. From preliminary 
runs |Q] we estimated this additive mass renormalization to be around 0.02. We ran the present 
simulations with lattice quark masses am uc i = 0.025, am s = 0.103. To determine the actual additive 
mass renormalization we calculated ratios 

P *®= 2(01^(^(0)10) ' = (5 " 1} 

For both correlator functions we used low-mode averaging and found well-defined plateaus, typ- 
ically beyond t « 5. These we denote by m^ 1 . Plotting m^fc a g a i nst m q = (m a + mi,)/2 and 
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Figure 4: Left: The chiral extrapolation of m AWI (m c/ ) using the combinations ud, us, and ss for the quark 
flavors. We find an additive mass renormalization of amo = 0.0196(2) and a renormalization constant of 
Z = 0.8837(7). Right: The low-lying eigenvalues of the massive Dirac operator for am u = 0.025. The red 
line shows the quark mass with this mass shift subtracted, am su b = am q — aniQ. 



extrapolating in m q to m AWI = we obtain the mass shift mo. Using naive currents, i.e. the non- 
conserved axial-vector current in the enumerator and the pseudoscalar density in the denominator 



of Eq. (5.1), we also get a multiplicative renormalization factor Z: 

1 



m^ 1 =Z(m q -m ) , m q = ^(m a + m h ) . (5.2) 



As seen on the left hand side of Fig. ^ the actual values are described very well by this linear 
dependence. We also plotted on the right hand side of Fig. || the 48 low-lying eigenvalues we stored 
for several of our configurations with am u = 0.025 on our 8 3 x 24 lattice. We find that they are lying 
approximatively on a Ginsparg- Wilson circle shifted by the subtracted mass am su b = am u — amQ = 
0.0054(2) (corresponding to 7.2(4)MeV) away from ReA = (red line). The subtracted lattice 
mass for the strange quark is am s — amo = 0.083, corresponding to m 1 lOMeV. 

5.2 Effective masses 

For the point-like pseudoscalar density P = xj/js W we calculated the zero-momentum (PP) 
correlator for three different quark flavor combinations (uu,us,ss) and found very stable effective 
mass plateaus shown in Fig. |[ 

Linear chiral extrapolation of m^ n in m q to the chiral limit, shown in Fig. |||, gives: 

(am eff ) 2 (m ) = 0.367(6) . (5.3) 



Using the finite volume mass gap in lowest order of chiral perturbation theory, Eq. (3.2), this 
leads to F = 92.7(4) MeV which agrees with the physical value of the pion decay constant, F n = 
92.4 MeV. However, this agreement is presumably accidental, since the spatial extend of our small 
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Figure 5: Mass plateau in the (Pab(0)Pba(t)) correlator, for ab = ud,us,ss. 
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Figure 6: The chiral extrapolations in the pseudoscalar channel. On the right hand side we extrapolate 
(oOTeff) 2 linearly in am q — + aniQ and on the left hand side we do the same for am e ff. 



lattices is uncomfortably small, and sizeable corrections to the leading order result of Eq. ( p.2| ) are 
expected. 



6. Summary and Outlook 

We calculated the autocorrelation for the 8 3 x 24 lattice and found autocorrelation times To ~ 
40 for plaquettes on smeared links and To ~ 80 for the (PP) correlation functions at fixed time 
distances. 

In the measurement of the meson correlators the low-mode averaging was found to improve 
the signal significantly. Since the low lying eigenmodes of the Dirac operator were available on 
each configuration from our updating, this did not require extensive computing resources. 

We determined the additive mass renormalization and verified that our light quark masses are 
near to the chiral limit. We found that the AWI mass depends nearly linearly on the average of the 



6 



First results in QCD with 2+1 light flavors using the fixed-point action 



Dieter Hierl 



corresponding lattice quark masses. We found that the eigenvalues lie in good approximation on 
the shifted Ginsparg-Wilson circle. 

For the pseudoscalar spectrum we obtained the finite volume mass gap in the chiral limit. 
Comparing it to the lowest order theoretical prediction we obtained a value for the low energy 
constant F. The surprising agreement with the physical value is presumably an accident since 
finite volume effects could be quite large at this lattice size. 

In the near future we will also have access to larger lattices, on which we plan to analyze 



conserved currents [ffl] and all-to-all propagators [11, 12] 
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